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Abstract. Wavelet analysis offers an alternative to Fourier based time-series 
analysis and is particularly useful when spectral features are time dependent. 
We analyze the solar neutrino capture rate of the radiochemical Homestake 
chlorine experiment with abbreviated Morlet wavelets, using Foster's (AJ, 
111,1709(1996)) rescaled wavelet technique. We emphasize the complemen- 
tarity of wavelet analysis to Fourier analysis. Wavelet analysis confirms the 
results of previously undertaken Fourier analysis. The Homestake data seem 
to contain a harmonic content with periodicities of 4.76 yr, 1.89 yr, 0.85 yr, 
and 0.51 yr. Wavelet analysis reveals that the 4.76 yr and 1.89 yr periods 
show an almost constant behavior over the 25 yr Homestake data record, 
while the 0.85 yr and 0.51 yr periods exhibit a transient phenomenon. The 
analysis does not show strong evidence for a period of the solar 11 yr cycle. 
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1. The Solar Neutrino Problem 

Solar neutrinos were detected by four experiments, the radiochemical 
Homestake chlorine experiment, the Kamiokande water Cherenkov experi- 
ment, and the two radiochemical gallium experiments, GALLEX and SAGE. 
In these four experiments, typically less than or of the order of 50 neutrino 
events were observed per year. These experiments confirmed that the sun 
shines by nuclear fusion reactions among light elements, burning hydrogen 
into helium, and that solar neutrinos have been observed in approximately 
the number and with the energies expected (Bahcall 1989, 1996). In April, 
1996, the Super-Kamiokande experiment began to operate and inaugurated 
a new area of high-statistics tests of the currently widely accepted standard 
solar models and standard electroweak theory (Suzuki 1997). 

The four experiments - Homestake, Kamiokande, GALLEX, and SAGE 
- have all observed solar neutrino fluxes with intensities that are within a 
factor of a few of those predicted by standard solar models. Three of the 
experiments (Homestake, GALLEX, and SAGE) are radiochemical and each 
radiochemical experiment measures one number, the total rate at which neu- 
trinos above a fixed energy threshold, which depends on the detector, are 
emitted by the sun. The only electronic (real-time, non-radiochemical) de- 
tector among the four experiments, Kamiokande, has shown that the neutri- 
nos come from the sun, by measuring the recoil directions of the electrons 
scattered by solar neutrinos (Bahcall 1996). 

Despite continual improvements of numerical standard solar model cal- 
culations of neutrino fluxes over more than three decades, the discrepancies 
between the different solar neutrino experiments and standard solar model- 
ing have increased with time. All four of the solar neutrino experiments yield 
neutrino capture rates that are significantly less than predicted by standard 
solar models. 

Solar modeling is required in order to predict the rate of nuclear fusion 
by the proton- proton chain reactions. In a standard solar model, about 
99% of the energy generation is produced by reactions in the proton-proton 
chain. The most important neutrino producing reactions are the low-energy 
pp (continuum: — OAMeV), pep (one line: lAMeV), 1 Be (two lines: 
0.86MeV, 0.38MeV), and the high-energy 8 B (continuum: - UMeV) 
neutrinos. 

In the simplest version of standard electroweak theory, neutrinos are 
massless and neutrino flavors are separately conserved; neutrinos do not os- 
cillate or decay to neutrinos with a different lepton number or energy. The 
minimal standard electroweak theory has been confirmed widely in precision 
laboratory tests (Barnett et al. 1996). Particle physics generalizations of the 
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standard electroweak theory (including left-right symmetry, grand unifica- 
tion, and supersymmetry) suggest that the neutrino may have a mass in the 
range 10~ 6 eU < m v < lO~ 2 eV. Precision laboratory particle tests have not 
yet substantiated such a suggestion. 

The first solar neutrino detection experiment to be performed was the 
crucial radiochemical Homestake chlorine experiment, which detects elec- 
tron neutrinos that are more energetic than 0.81MeV. After almost three 
decades of operation of this experiment, the measured neutrino capture rate 
is 2.56 ± 0.22 SNU, which is a factor ~ 3.6 less than is predicted by the most 
detailed theoretical calculations, 9.3jz};l SNU (a SNU is a convenient unit 
to describe the measured rates of solar neutrino experiments: 10~ 36 interac- 
tions per target atom per second). The predicted neutrino capture rate in the 
Homestake experiment is dominated by the rare, high-energy 8 B neutrinos, 
although 7 Be neutrinos also contribute significantly. According to standard 
model calculations, the pep neutrinos and the CNO neutrinos are expected 
to contribute less than 1SNU to the total neutrino capture rate. The discrep- 
ancy between the standard solar model calculations and the observations of 
the radiochemical Homestake chlorine experiment came to be known as the 
solar neutrino problem (Davis 1968, 1992, 1996; Zimmerman 1996). 

The comparison of the experimental results of the four solar neutrino 
experiments shows (i) the smaller than predicted absolute event rates in 
the Homestake and Kamiokande experiments, (ii) the incompatibility of the 
Homestake and Kamiokande experiments, and (iii) the very low neutrino cap- 
ture rates in the GALLEX and SAGE experiments which may imply a great 
reduction of 7 Be neutrinos, although 8 B neutrinos are observed, with respect 
to the value predicted by the standard solar models. The conclusion based 
on these three facts is that (i) either at least three of the four solar neutrino 
experiments (GALLEX and SAGE and either Homestake and Kamiokande) 
have yielded misleading results, (ii) physics beyond the standard electroweak 
theory is required to change the neutrino energy spectrum or flavor content 
after the neutrinos are produced in the central region of the sun, or (iii) the 
solar neutrino flux is varying over time due to unknown physical phenom- 
ena. This paper is addressing the search for such a variation of the neutrino 
capture rate in the Homestake experiment. 

In Section 2 we review the current status of the search for periodicities in 
the Homestake data record and summarize the results obtained by Fourier 
and Lomb-Scargle analysis of the data. 

The relationship between Fourier analysis and wavelet analysis is presented 
in Section 3.1. 
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The method of the weighted wavelet Z-transform and weighted wavelet am- 
plitude is outlined in Section 3.2. This special implementation of wavelet 
analysis takes into account the fact that the Homestake data are unevenly 
spaced, contain spikes and gaps in the record, and cover only a finite interval 
of time spanning 25 years. 

The results for the search for periodicities based on wavelet analysis are 
given in Section 3.3. Wavelet analysis confirms that the Homestake data 
contain a series of periodicities of 4.76 yr, 1.89 yr, 0.85 yr, and 0.51 yr. The 
time dependence of these spectral features reveals that the first two periods 
exhibit an almost constant behavior over the 25 years of Homestake data 
record, while the latter two periods show a transient nature with variability 
of the actual value of the period. We find no strong evidence for a period of 
the order of 11 yr in the data. 

Conclusions are discussed briefly in Section 4. 

2. Is the Solar Neutrino Flux Constant Over Time ? 

The solar neutrino flux has been inferred from the neutrino capture rate in the 
Homestake experiment (Figure 1), measured over the past 25 years (1970.281- 
1994.388) in 108 runs (Davis 1968, 1992, 1996). Recently, the Homestake 
data have been completely reanalyzed by the operators of the experiment, 
leading to changes in the uncertainties of the measurements (Davis 1996). 
We emphasize this fact by dubbing them "new" solar neutrino capture rate 
data. The observed neutrino capture rate is, on average, several times smaller 
than the predicted one based on standard solar models. In the following we 
are not concerned with the discrepancy between the average rate and the 
predicted rate, but with a possible time dependence of the neutrino capture 
rate in the radiochemical Homestake chlorine experiment. Suggestions that 
the neutrino capture rate is anticorrelated with solar activity, reflected by 
the number of sunspots, were made by Davis (1992, 1996), in particular 
when the neutrino capture rate went through an apparent minimum around 
1980 (coinciding with the time when sunspot cycle 21 reached its maximum). 
Such a suggestion was particularly encouraged when analyzing the five-month 
moving average of the Homestake data record. 

Bahcall, Field, and Press (1987) and Bahcall and Press (1991) showed 
that conclusions about correlation with sunspot number, essentially based 
on the Spearman rank order correlation coefficient, Kendall's tau, and the 
"shuffle test", change significantly when one uses either the upper errors or 
the average errors of individual runs. In the following we have chosen to limit 
the analysis of the solar neutrino capture rate in the radiochemical Homestake 
chlorine experiment to the average value for each run, although the analysis of 
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Bahcall, Field, and Press (1987) and Bahcall and Press (1991) indicated that 
the size of the error bars for the runs with a very small neutrino capture rate 
is important for determining the time dependence of the solar neutrino flux. 
Sturrock, Walther, and Wheatland (1997) presented statistical tests for the 
variation of the neutrino capture rate in the Homestake experiment and find 
evidence that the rate is not constant over time. Particularly, they searched 
for a low-frequency content in the Homestake data and concluded, through 
spectrum analysis based on the use of maximum likelihood estimation, that 
there is no evidence for periods of the order of 11 yr (solar cycle) but some 
evidence for periods of 2.1 yr (Sakurai periodicity; Sakurai 1979) and 0.44 
yr (Rieger periodicity; Rieger et al. 1984). We do emphasize here that at 
the current point of time the latter two periodicities can not be related to 
or explained by any solar periodic process in the central region of the sun 
based on standard solar modeling (Sakurai presented evidence that such a 
quasi-biennial periodicity is contained in the Homestake data; Rieger et al. 
have shown that the occurrence rate of solar flares exhibits a periodicity of 
about this length). Sturrock, Walther, and Wheatland (1997) also searched 
the Homestake data for high-frequency content, i.e. variations of the order of 
or less than the average exposure time of individual runs, and find support 
for an earlier proposition (Sturrock and Bai 1992) that the neutrino capture 
rate may be modulated at a frequency that could be related to the sidereal 
rotational frequency of 13.88 cycles yr -1 of the sun's radiative zone. 

Standard solar modeling predicts that the neutrino fluxes are constant in 
time. Neither solar physics nor nuclear physics of the standard solar model 
allow variation of the solar neutrino flux over time. Additionally, the four 
solar neutrino experiments are low event rates experiments with large error 
bars on each of the measured values. Despite of this, it can not be excluded, a 
priori, that unknown physical phenomena may lead to a variation of the flux 
additional to seasonal variations caused by the Earth's orbital eccentricity. 

The simplest technique available for investigating periodicities in a data 
record is the Fourier analysis; i.e., the comparison between a record and a si- 
nusoidal signal with a given frequency. Fourier analysis has a long tradition 
in analyzing "quantities subject to irregular fluctuations" (Einstein 1914), 
however, the traditional Fourier transform method often may not be suitable 
for representing the real structure of sparse data (Foster 1996a; Sturrock, 
Walther, and Wheatland 1997). This shortcoming is evident when suspected 
periods and their amplitudes change with time and the sampling extends over 
a finite interval. The application of Fourier analysis to such a record may lead 
to the identification of spurious periodicities. Therefore, a method of avoiding 
the identification of artificial periodicities is needed, and wavelet transform 
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is a technique that seems to suite this purpose (Strang 1993; Jawerth and 
Sweldens 1994). The wavelet method can reveal variations of spectral fea- 
tures in a time dependent solar neutrino flux. Foster (1996b) introduced a 
rescaled wavelet analysis technique that appears to be an effective method of 
reducing artificial periodicities which may result from Fourier analysis. Three 
characteristics of the Homestake record complicate period analysis. The first 
is the limited duration of the record, the second is the presence of spikes 
and gaps in the record, and the third is the uneven distribution of the data. 
Gaps are related to physical and practical constraints; i.e., the mean-half life 
of argon in the radiochemical Homestake chlorine experiment ( 35days) and 
facility maintenance. Standard wavelet methods require time series to be 
regularly distributed in time. Extended gaps in the Homestake record can 
not be filled by interpolation because of the unknown pattern of the signal 
itself. Commonly, these constraints appear in the context of both Fourier 
and wavelet analysis. An extension of Fourier technique that attempts to 
circumvent one of those limitations is the Lomb-Scargle periodogram ( Lomb 
1976; Scargle 1982). Its aim is to correct the trial functions (cosut, smut) of 
the Fourier transform to preserve its normalization condition on an uneven 
grid of times of a data record. For this purpose, an extraction of the mean 
value of the signal under investigation and a phase shift of trial functions are 
used. The appealing feature of the Lomb-Scargle normalized periodogram is 
that it weights data on a per-point rather than a per-time-interval basis. Ta- 
ble 1 summarizes the results of previous Fourier and Lomb-Scargle analysis 
of the Homestake data record, where the uneven distribution of the data had 
been taken into account. 

In a way similar to that in which the Lomb-Scargle periodogram comple- 
ments the Fourier analysis, the wavelet analysis has been extended by Foster 
(1996b) to the rescaled wavelet technique. The rescaled wavelet technique 
can be applied to short data series with spikes and gaps. The idea behind 
this technique is basically the same as that behind the Lomb-Scargle peri- 
odogram: wavelet trial functions are extended and corrected by a rescale 
procedure to satisfy the admissibility condition on the uneven grid of times 
of a data record. In this technique the wavelet transform is understood as 
a projection onto trial functions which do resemble the shape of the data 
record. 

We describe the rescaled wavelet technique and then compare the results 
obtained from the rescaled wavelet transform technique to previous period- 
icity searches performed on the Homestake record with Fourier analysis and 
the Lomb-Scargle normalized periodogram method (Table 1). Wavelet anal- 
ysis allows the study of the variability of spectral features over time unlike 
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the Fourier transform. 



3. Wavelet Analysis of the Homestake Data Record 

3.1. Wavelet and Fourier Spectra 

Fourier analysis fails when one needs to consider the time dependence 
of spectral features. By contrast, wavelets are efficient in multiscale analysis 
(Strang 1993; Jawerth and Sweldens 1994). They have a localized, oscillating 
form so that, unlike sinusoids, they are localized near time r and decay if 
\t — r\ exceeds a characteristic scale, a . Therefore, the wavelet representation 
can be considered to be a mathematical microscope with variable position 
and magnification. The wavelet transform represents one-dimensional signals 
as a function of both time and frequency (position and scale) and is similar 
to a local, filtered Fourier transform that can be obtained by dilating and 
translating the wavelet and then convolving it with the signal (Walker 1997). 
We have a real signal / of a real variable t , represented by the data record 
f(t). One defines its Fourier transform F{k) and its Fourier spectrum Pf(&) 
accordingly, 

/+oo 
dtf(t)e~ m , (1) 
-oo 

P F (k) = ^\F(k)\ 2 for k>0. (2) 
The total energy E of the signal / is defined such that 

1 r+oo 1 r+oo r+oo 

E = 2 loo dm)l = 4^ loo dkimi = Jo dkMk) - (3) 

The function ijj is called an analyzing wavelet, if it verifies the admissi- 
bility condition 

r+oo 

c^= dcou" 1 ^^)] 2 < +oo, (4) 
Jo 

where 

$(u) = dti){t)e- ibJt (5) 

J — oo 

is the Fourier transform of the wavelet. Condition (4) implies that the in- 
tegrand defining should be integrable at oj = and therefore that the 
wavelet has a zero mean, "0(0) = 0. The wavelet will also decay to zero as 
the variable tends to zero. A stronger condition is to require cancellation up 
to some order p, that is 

/+oo r+oo 
dtt n ifj{t) = for n = 0,1,..., p- 1, and / dtt p if){t) ^ 0. (6) 
-oo J —oo 
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If the function ijj is integrable and square integrable, then the admissibility 
condition (4) is equivalent to a cancellation condition of the order of at least 
zero and condition (6) is equivalent, up to a multiplicative constant, to the 
existence of a bounded continuous function <fi with ^(0) = 1 and 0(oo) = 
such that 

${u>)= u?4>{u). (7) 

Additionally, the wavelet ip should be localized both in physical and in Fourier 
space (time and frequency). This requirement means that the time spread, 
At, and the frequency spread, Acu, of ip must satisfy the Heisenberg relation 
AtAcu = const., where the constant cannot be smaller than 2ty. One defines 
the wavelet transform W(a, r) of the signal /, 

1 r+oo t — t 
W{a,r) = ^ m dtf(t)P( ), (8) 

CI ' J —oo CI 

where a and r denote the dilation (frequency) and translation (position) scal- 
ing factors, respectively. For \a\ << 1, the wavelet is a highly concentrated 
shrunken version of ip with frequency content mostly in the high frequency 
range. If \a\ >> 1, the wavelet is spread out and contains mainly low fre- 
quencies. In practice, the integral limits may be replaced by values at which 
the wavelet amplitude is below some threshold, and t is evaluated at data 
points of the time series, which need not be the same times as r . The con- 
tinuous wavelet transform can also be computed from the Fourier transform 
of the signal /: 

(7 1 / 2 r+oo 

W(a,r) = / dcoF(co)^*(au)e tuJT . (9) 

2lT J-oo 

From the energy conservation property of the wavelet transform, 

1 r+oo An r+oo 

we define the local wavelet spectrum 

2 

for k > 0, (11) 

where k$ denotes the peak frequency of the analyzing wavelet -0. The local 
wavelet spectrum measures the contribution to the total energy coming from 
the vicinity of the point t and frequency k, where the vicinity depends on 
the shape of the analyzing wavelet in physical and Fourier space. From the 
local wavelet spectrum we derive a mean wavelet spectrum Pw{k): 

/+oo 
dtP w (k,t), (12) 
-oo 



Pw(k,t) 



2c^k 
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which is directly related to the total energy E of the signal / 

r+oo 

E= / dkP w (k). (13) 
Jo 

Equations (2), (9), (11), and (12) reveal the relationship between the 
Fourier spectrum Pf(uj) and the mean wavelet spectrum Pw(k), namely, 

(14) 

Equation (14) shows that the mean wavelet spectrum is the average of the 
Fourier spectrum weighted by the square of the Fourier transform of the 
analyzing wavelet tp shifted at frequency k. The higher the frequency k 
the wider the averaging interval. Based on this behavior of the mean wavelet 
spectrum it can be shown that the mean wavelet spectrum at high frequencies 
depends on the behavior of the analyzing wavelet at low frequencies. 

In the following we use the family of analyzing wavelets which are si- 
nusoids in a Gaussian envelope, known as modified Morlet wavelets. This 
family is well suited for adjusting the tradeoff between time and frequency 
resolution, 

ip(z, n) = (e iz - e'^^j (15) 

2 2 

where z = — and e ~ is a small correction term which has been 

a 

inserted (modifying the Morlet wavelet) so that it satisfies the admissibility 
condition (4). The parameter in the corrective term is related to Foster's 
(1996b) constant c by 2n 2 Tx 2 = c _1 , which suits better the wavelet used in 
this analysis; it determines how rapidly the analyzing wavelet decays. For 
n ~ 1, the parameter n corresponds to the number of sinusoidal periods 
which fit between the inflection points of the Gaussian envelope. For n ~ 2, 
the corrective term can even be neglected as it becomes very small. The 
larger the value of n the wider the window and the better the resolution in 
frequency but the worse the resolution in time. In the following numerical 
wavelet analysis we will adopt the value n = 2 (c ~ 0.0127), for this choice 
the corrective term is negligible and leads to the abbreviated Morlet wavelet 

ip(z,n)=e lz ~^^=e t ~ £w . (16) 

The chosen value n = 2 fine-tunes the decay rate of the wavelet. 

Employing the abbreviated Morlet wavelet leads to the abbreviated Mor- 
let transform and, for the discrete data record under consideration, it defines 
the discrete wavelet transform. The abbreviated Morlet transform looks sim- 

(t-r) 2 

ilar to a windowed Fourier transform, with the window e 2^7^ ; exhibiting 



Pw(k) 



c^k Jo 



dujPpiuj) 



i( koUJ \ 
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the factor a~ 2 in the exponential that makes the size of the window fre- 
quency dependent. Assuming, for example, the total width of the abbrevi- 
ated wavelet is about 7 years, we can find the correlation between this curve 
and the first 7 years of the Homestake data record shown in Figure 1. This 
single number gives a measure of the projection of this wavelet on the data 
record during the 1970 - 1977 period; i.e. how much (amplitude) does our 7 
year period resemble a sine wave of this width (frequency). By sliding this 
wavelet along the Homestake data record one can then construct a new data 
record of the projection amplitude versus time. Eventually, one can vary 
the scale of the wavelet by changing its width. This is the real advantage of 
wavelet analysis over a moving Fourier spectrum. For a window of a certain 
width, the sliding Fourier transform is fitting different numbers of waves; i.e. 
there can be many high-frequency waves within a window, while the same 
window can only contain a few (or even less than one) low-frequency waves. 
The wavelet analysis uses a wavelet of the exact same shape, only the size 
scales up and down with the size of the window. 

3.2. Weighted Wavelet Z-Transform and Weighted Wavelet Amplitude 
The discrete analogue to (8) is 

1 N /f — t\ 

for the observed time series consisting of N data values f(t a ), taken at a 
discrete set of N times t a , a = 1, 2, N. For the abbreviated Morlet wavelet 
(16) the real and imaginary parts of the discrete wavelet transform are 

Re(W) = ^72 E f(t a )e-£^cos ( ^— -) (18) 

and 

Im{W) = f] f(t a )e-^sin -) . (19) 

When treating the discrete wavelet transform (17) as a projection, it can be 
interpreted as a weighted projection onto the trial function 

$(t) = e^-ir (20) 

with the statistical weights chosen as 

w a — e 2» 2 " 2 » 2 . (21) 
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In order to tackle the problem of uneven time spacing with the discrete 
wavelet transform, similar to a procedure used in discrete Fourier tansform, 
Foster (1996b) suggested to include a third trial function, the constant func- 
tion l(t) = 1 for all t. Therefore, a weighted projection onto three trial 
functions is performed, namely 

*i(t) = Ht), (22) 
$ 2 (t) = cos(— ), (23) 

(X 

$ 3 (t) = s i n (— ). (24) 
ct 

Projection determines the coefficient y a of a set of r trial functions $ a (t), a = 
1, 2, . . . , r for which the function under consideration, 

l/(*)=E !/■(*)*«(*). ( 25 ) 

a 

fits best the data record. Such a projection can be computed by defining the 
inner product of two functions u(t) and v(t) as 

<u\v>= ^ W f^ )v ^\ (26) 

£/J=l w /3 

where w a denotes the statistical weight assigned to the data point a. The 
best-fit coefficients of the trial functions are determined by multiplying the 
inverse of the S-matrix, that is the matrix of inner products (26) of the trial 
functions (S a b =< & a \&b >) by the vector of inner products of the trial 
functions with the data record 

y. = E^<^l/>- ( 27 ) 
b 

Based on the above procedure, Foster (1996a) defines the power for evaluating 
the projection statistically as 



N 



l)s 2 



< >< <W > - < 1|/ > 2 ) (28) 



for the set of trial functions, where N is the number and s 2 the estimated 
variance of the data record. The power (28) is a chi-square statistic with 
r — 1 degrees of freedom and expected value 1. 

In the case of a weighted projection, taking into account (20) and (21), 
the number of data N in (28) has to be replaced by an effective number 
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(Foster 1996a, 1996b) 



N, 



eff 



(E<) 



(29) 



and for the variance s 2 in (28) the weighted estimated variance 

2 N eff V f 



N, 



eff 



(30) 



has to be used. Accordingly, in (30) Vf denotes the weighted variation of the 
data record, 



Vf 



Eg W a P(t a ) 



Eg Wafjtg) 

Ea^a 



-<f\f> ~<l\f>' 



(31) 



and V y is the weighted variation of the function y 



Vy = 



Eg w a y 2 (t c 
Ea^a 



Eg w a y(t c 



1 2 



=< y\y > - < l\y > 2 



Ea^a 

Thus, Foster's (1996b) weighted wavelet transform (WWT) is defined by 

(Neff - l)V y 



(32) 



WWT = 



2V f 



(33) 



for fixed a and r, as a chi-square statistic with two degrees of freedom and 
expected value one. 

The weighted wavelet transform in equation (33) depends sensitively on 
the effective number of data N e ff which leads to a shift of the WWT peaks 
to lower frequencies because at lower frequencies the window is wider and 
thus more data points are sampled and the effective number N e ff is larger. 
To compensate for this fact, a test statistics which is less sensitive to the 
effective number of data can be applied. Such a statistics for projections 
was suggested by Foster (1996a), which leads to the weighted wavelet Z- 
transform (WWZ) 

wwzJ -w^w (34) 

The statistical behavior of (34) is derived for a projection where the statistical 
weights are inversely proportional to the variance of the data. 

Although the WWZ may be an excellent locator of the signal frequency, 
it is a poor measure of amplitude. With a projection it is easy to define 
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the amplitude of the corresponding periodic fluctuation: it is the square 
root of the sum of the squares of the expansion coefficients for the sine and 
cosine functions, which Foster (1996b) calls the weighted wavelet amplitude 
(WWA), 

WWA={{y 2 f + {y,f) l '\ (35) 

Therefore the weighted wavelet transform allows us to compute the amplitude 
by the WWA, after determining the period from the WWZ. 

3.3. Results of the Wavelet Analysis 

For the Homestake data we have subtracted the average value from the 
data prior to analysis. The WWZ was computed over the frequency range 
0-2.5 cycles yr -1 in steps of 0.01 cycles yr -1 . It must be noted that we 
do read off the significance of a peak directly from the heights of the peaks 
in the WWZ. Fig. 2 shows the perspective wavelet map for the analysed 
Homestake data. 

In examining the WWZ, as seen in Fig. 3, in the neighborhood of the 
frequency of the solar cycle (0.09 cycles yr" 1 ) we find no strong evidence of a 
peak, even if the WWZ shows a peak at 0.12 cycles yr -1 . This is compatible 
with the analysis of Haubold and Gerth (1990), Haubold and Beer (1991), 
and Sturrock, Walther, and Wheatland (1997). The examination of the 
spectrum in the neighborhood of the Sakurai quasi-biennial periodicity (0.45 

- 0.50 cycles yr -1 ) and the Rieger periodicity (supposed to be located at 2.30 

- 2.40 cycles yr -1 ) we find there are peaks with high power at frequencies 0.53 
and 1.93 cycles yr -1 , respectively. However, the power spectrum reveals even 
stronger peaks at frequencies 0.21 cycles yr -1 and 1.18 cycles yr -1 which can 
not be related to any known solar periodic process in the central region of 
the sun at this point of time. The most prominent peaks in the WWZ are 
the frequencies at 0.21 and 0.53 cycles yr" 1 . Inspection of the development 
of the WWZ in Fig. 4 over the 25 yr of Homestake record indicates that 
the value of the frequencies at 0.21 cycles yr -1 and 0.53 cycles yr -1 are 
stable at their positions in the WWZ, showing, however, strong variation 
in the value of the WWZ. Contrary, the frequencies 1.18 cycles yr -1 and 
1.93 cycles yr -1 exhibit a wide variation of their values in the WWZ over 
25 yr of Homestake record, pointing towards a transient phenomenon. From 
inspection of Figures 5 and 6, showing the Fourier spectrum and amplitude of 
the Homestake data, respectively, computed over the same frequency range 
with the same frequency step as the WWZ, we find evidence for peaks near 
0.21 cycles yr -1 and 0.53 cycles yr -1 as well as near 1.18 cycles yr -1 and 1.93 
cycles yr -1 which have been also revealed in the wavelet analysis. 

Examination of the WWZ and power spectra shows that a series of peri- 
odicities (0.21, 0.53, 1.18, 1.93 cycles yr" 1 ) shows up in each analysis. 



13 



4. Conclusion 

In a way similar to that in which Lomb-Scargle periodogram complements 
Fourier analysis, Foster's weighted wavelet Z-transform extends wavelet anal- 
ysis. In addition, the WWZ allows to study the variability of periodicities 
over time, unlike Fourier analysis. The computation of the WWZ confirms, 
to a reasonable extend, the results of Lomb-Scargle and Fourier analysis, 
summarized in Table 1, applied to the unevenly spaced data yielded by the 
radiochemical Homestake chlorine experiment, in that the data have a har- 
monic content. Analysis, using the abbreviated Morlet wavelet, indicates 
that this harmonic content is represented by a series of periodicities of 0.21, 
0.53, 1.18, and 1.93 cycles yr -1 . 

The two "low-frequency" values of them (0.21 and 0.53 cycles yr -1 ) show a 
stable behavior over the analyzed 25 years of Homestake record while the 
latter two frequencies exhibit a transient behavior by large variation of their 
values over the 25 yr of data analyzed. Two periodicities correspond with 
Sakurai's periodicity (0.45 - 0.50 cycles yr -1 ) and Rieger's periodicity (2.30 
-2.40 cycles yr -1 ). None of the above analysis shows strong evidence for a 
variation in the Homestake data with the 11 yr solar cycle. 

A variation of the solar neutrino capture rate in the Homestake experi- 
ment may result from fluctuations of nuclear fusion or by a modulation of 
the capture rate due to the solar magnetic field or helioseismic waves but 
both phenomena are not supported by standard solar modeling (Sturrock, 
Walther, and Wheatland 1997). 
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Table and Figure Captions 

Table 1. Summary of periods ir revealed by Fourier, Lomb-Scargle, and 
wavelet analysis of the Homestake record as discussed in the text and in 
Haubold and Gerth (1990) and Haubold and Beer (1991). 

Fig. 1. A plot of the individual solar neutrino capture rates from the Home- 
stake experiment, runs nos. 18 to 133 (1970.281 - 1994.388) (Davis 1996). 

Fig. 2. The perspective wavelet map (frequency-time- WWZ) for the Home- 
stake record. The value of WWZ is approximately an F-statistic with N e ^ 
(the effective number of data for the given time and frequency being ana- 
lyzed; eq (29)), two degrees of freedom, and expected value one. The graph 
indicates whether or not there is a periodicity at a given time of the given 
frequency. 

Fig. 3. The value of WWZ from Fig. 2 as a function of frequency showing 
the series of periodicities in the Homestake record as discussed in the text. 

Fig. 4. The value of the WWZ from Fig. 2 showing the time evolution of 
the amplitude of the periodicities in the Homestake record. 

Fig. 5. Spectrum of the Homestake data for low-frequencies, obtained by 
Fourier analysis. 

Fig. 6. Amplitude of the Homestake data, as determined by Fourier analysis, 
corresponding to the spectrum shown in Fig. 5. 





n[yr] 

No. runs 


71 < 1 


1 < 7T < 2 


2 < 7T < 3 


3 < 7T < 4 


Fourier 
18-69 

(1970-1981) 


0.7 
0.5 


1.63 
1.30 


2.14 


3.00 


Fourier 
18-89 

(1970-1985) 


0.83 
0.61 
0.54 
0.51 


1.61 


2.13 




Lomb-Scargle 

18-109 

(1970-1990) 


0.7 
0.54 


1.30 






Fourier 

18-133 

(1970-1994) 


0.7 
0.55 
0.53 
0.51 


1.75 
1.59 


2.04 




Wavelet 

18-133 

(1970-1994) 


0.85 
0.51 


1.89 






n[yr] 

No. runs 


4 < 7i < 5 


5 < 7T < 6 


8 < 7T < 9 


9 < 7T < 10 


Fourier 
18-69 

(1970-1981) 


4.90 




8.33 




Fourier 
18-89 

(1970-1985) 






8.33 




Lomb-Scargle 

18-109 

(1970-1990) 


4.80 






9.60 


Fourier 

18-133 

(1970-1994) 


4.55 








Wavelet 

18-133 

(1970-1994) 


4.76 
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